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Abstract 

Modelling the chaotic states in terms of the Gaussian Orthogonal Ensemble of 
random matrices (GOE), we investigate the interaction of the GOE with regular 
bound states. The eigenvalues of the latter may or may not be embedded in the 
GOE spectrum. We derive a generalized form of the Pastur equation for the average 
Green's function. We use that equation to study the average and the variance of the 
shift of the regular states, their spreading width, and the deformation of the GOE 
spectrum non-perturbatively. We compare our results with various perturbative 
approaches. 

Key words: Chaos, Random-Matrix Theory, Doorway State, Spreading Width 
PACS: 21.10.-k, 21.60.-n 



1 Introduction 



Predictions of the Gaussian Orthogonal Ensemble of Random Matrices (GOE) 
offer the best first guess of spectral fiuctuation properties of a system about 
which nothing is known beyond the fact that it is time-reversal invariant [Ij. 
Such a description is completely adequate if the system is classically chaotic. 
Sometimes, additional but incomplete dynamical information exists which 
calls for an extension of the GOE. To make the need for such an extension 
plausible, we mention several examples, (i) In spherical atomic nuclei, the shell 
model is known to provide a proper dynamical description of the low-lying 
states. Near neutron threshold, however, the density of states in medium- 
mass and heavy nuclei is so high and the levels are of such complexity that 
the shell-model approach becomes unfeasible (even though it is probably still 
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adequate in principle). Moreover, there is evidence that the spectral fluctua- 
tion properties of these states coincide with GOE predictions. The states are 
accordingly modelled in terms of the GOE. The various realizations of the 
GOE may then be viewed as corresponding to different choices of the residual 
interaction of the shell model. The proper (i.e., realistic) choice is unknown. 
How do such different choices affect the predictions of the shell model at lower 
excitation energies? The question can be explored by studying an extension 
of the GOE comprising in addition to the GOE Hamiltonian a set of discrete 
states with energies below the continuous GOE spectrum. The states interact 
with each other and with the GOE Hamiltonian. (ii) In condensed-matter the- 
ory, a quantum dot carrying a single state may be coupled to a finite reservoir 
of interacting electrons. The reservoir may be modelled by a GOE Hamilto- 
nian. The single state on the dot lies within the GOE spectrum. How is the 
state affected by the interaction with the GOE? That same question arises 
in finite Fermi systems when a particular mode of excitation, which is not 
an eigenstate of the Hamiltonian, occurs in a sea of states carrying the same 
quantum numbers. Lacking a better description, those states are modelled in 
terms of the GOE. In nuclear physics and in quantum chemistry, the single 
state is referred to as a doorway state. 

To address these and related questions, we study in the present paper the 
spectral properties of a generic system consisting of a GOE Hamiltonian in- 
teracting with a set of discrete states. The states may or may not lie within 
the GOE spectrum. Related problems have been studied by several authors. 
The doorway state mechanism has found wide interest. Early reviews may be 
found in Refs. [2|3] . see also Ref. The interaction between the higher-lying 
complex states and the low-lying states of the shell model was, for instance, 
addressed in a series of papers by Feshbach and collaborators |5] and taken 
up again in Refs. [B]. The present work goes beyond these papers in that we 
present a comprehensive treatment of the problem within the framework of 
random-matrix theory. 

The paper is organized as follows. The random-matrix model is defined in 
Section [2J A straight perturbative approach which is useful for purposes of 
orientation is described in Section [31 both for non-overlapping and for over- 
lapping spectra. Section H] is the central theoretical piece of the paper. We 
derive and analyze an extension of the Pastur equation for the average Green's 
function of the system in the limit of infinite matrix dimension of the GOE 
Hamiltonian. In Section O we specialize that equation to the case of a sin- 
gle state interacting with the GOE. We solve the generalized Pastur equation 
perturbatively. In Section [6] we present numerical results based upon an exact 
solution of the Pastur equation. These are compared with the various analyt- 
ical approximations derived earlier. Section [7] contains our conclusions. Some 
technical details are deferred to the Appendix. 
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2 Formulation of the Problem 



We consider a decomposition of Hilbert space into two subspaces defined by 
orthogonal projection operators P = P'^ and Q = Q\ with P^ = P, = Q, 
PQ = 0. The dimension of P-space is M, that of Q-space is A^. We shall as- 
sume that ^ 1 but (for reasons given below) will restrict ourselves to small 
values of M. The dynamics in P-space is determined by the Hamiltonian Hq. 
The motion in Q-space is assumed to be irregular or stochastic and described 
in terms of random-matrix theory. The Hamiltonian if^iOE _ qjjGOEq ^ 
member of the ensemble of Gaussian orthogonal random matrices (GOE). The 
matrix elements H^^^ with n,^ = 1, . . . , N are Gaussian-distributed random 
variables with zero mean value and second moment 

Here and in what follows, the average over the ensemble is denoted by angular 
brackets, and 2A denotes the radius of Wigner's semicircle, see Eq. Q) below. 
The Hamiltonians in P-space and in Q-space are coupled by an interaction 
V = PVQ + QVP. The total Hamiltonian has the form 

H = Ho + V + H^^^ . (2) 

We study how the stochastic dynamics in Q-space affects the dynamics in 
P-space. 



3 Perturbative Approach 



It is instructive to use first a perturbative approach as this gives some insight 
into the behavior of the system. We consider the case M = 1. The single 
state in P-space carries the index while the states in Q-space are labelled 
H = 1, . . . ,N. The matrix representation of the Hamiltonian ([2]) reads 



H 



(3) 



where Eq = {Hq)qq, and where Vq^ = Vf^o = are the matrix elements 
coupling the P-space to the Q-space. It is convenient to diagonalize H'^'^^. We 
call the eigenvalues i?^ and the transformed coupling matrix elements V^. The 
eigenvalues obey Wigner-Dyson statistics, and the V^'s are uncorrelated 
Gaussian random variables with zero mean value and with a common variance 
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V^. Moreover, the V^'s and the -E^'s are uncorrelated. Eq. ([3]) takes the form 

H=(- - 1. (4) 



3.1 Non-overlapping Spectra 



By definition, the GOE spectrum is centered a.t E = 0. We assume that the 
spectra of Hq and of H'^'^^ do not overlap. Hence, the distance |£^o| of the 
P-space state from the centre of the GOE spectrum must be larger than 2A, 
the radius of the semicircle, see Eq. below. The secular equation for H is 
easily found to read 

TV Y2 

We assume that (V^)^/^ is small in comparison with the difference between Eq 
and the closest end point of the semicircle. We accordingly write 

a = Eq + Sa (6) 

and assume that we may solve Eq. ([5]) by expanding in powers of 6a. This 
gives 

TV oo / 1 \ "+1 

-'^« = Ev^E^^ i^^r (7) 

f,=l n=0\EfM-EoJ 

which shows that Sa is of order V^. Keeping terms up to fourth order in 
we get 

AT T/2 / AT T/2 \ 



u 



We calculate both, the ensemble average {5a) and the variance {{5a — (5a;))^) 
of 5a. Different realizations of the GOE-Hamiltonian give rise to different 
values of 5a. The ensemble average of 5a yields the mean shift of Eq due to 
the interaction with the states in Q-space, and the variance of 5a is a measure 
of the fluctuation of the position of the state in P-space due to different 
realizations of the GOE Hamiltonian. The calculation of both quantities uses 
N ^ 1 and is sketched in the Appendix. We recall that for < 2 A, 

N 



p{E) = —^l-{El{2\)f (9) 

is the average density of states of the GOE, and that the usual deflnition of 
the spreading width for a state mixed with the GOE and located at P = 

ism 

= 27rVV(0) . (10) 
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With this definition, we find 



{Sa) 



27r 



/(xo) 



d/(xo) 
47rA dxo 



'111 



Here Xq = Eq/2\ obeys |xo| > 1, and 



I{xo) 



dx- 

1 X — Xo 



2Jx'o - 1 



-nxo 



arctan — , — arctan 



Xq 1 _ 



fl2) 



The second term on the right-hand side of Eq. (fTT!) suggests that the pertur- 
bation expansion proceeds in powers of T^/X. For |xo| ^ 1 we have /(xq) ~ 
— tt/ (2xo) and 

(13) 

2-C/o 

The shift is always away from the center of the semicircle and of order TK 
For the variance of Sa we obtain 



var(5a) 



iV \ 2^/ 



dxo 



I{xo) 



I{xo) 



d^ 



SttA 



dxo 



(14) 



Eq. f[T^ confirms the impression that the perturbation expansion proceeds in 
powers of T^/ X. Moreover we see that (var(5a))^/V (Sa) is of order A^-^/^^ ^^his 
is due to the fact that in Wick-contracting the V^'s, we reduce the number of 
independent summations over the Q-space states by one. For |xo| S> 1 and to 
lowest order in T^/X, we find 



vai{6a) 



1 rU 



VN2\Eo\ 

Eqs. ( [T3|) and (flSjl are in accord with the findings of Ref. [6]. 



(15) 



The interaction between the levels in P-space and those in Q-space induces 
an interaction amongst the levels in P-space. To estimate that interaction we 
consider two degenerate P-space states at energy Eq and use a representation 
of the GOE Hamiltonian in the form of Eqs. ([3]) and (jll). Perturbatively, the 
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matrix element of the induced interaction between these two states (labelled 
a and 6, respectively) is given by 

Using the same steps and notation as before, we find for the ensemble aver- 
age of V^'^'^ the approximate value {Va\Vb) / Eq. With (p the angle between the 
vectors Va^ and this is approximately equal to cos(0)(A/-Eo)(rV2)- Com- 
paring this expression with the average shift f|T3|) we see that the averaged 
induced interaction matrix element is smaller by the factor coscf). Because of 
the complexity of the states in Q-space, it is reasonable to expect that for 
^ 1 the two vectors Va and are approximately orthogonal. Then the 
induced interaction between two degenerate states in P-space is very small. If 
cos = the variance of V^'^ is small of order This case applies whenever 
we assume that the two P-space states are coupled with equal strength to the 
states in Q-space, see Sections [5] and El We conclude that it it reasonable to 
assume that the induced interaction between the states in P-space is small 
in comparison to the shift of each state. This is why we confine ourselves in 
Section [6] mainly to a one-dimensional P-space. 



3.2 Overlapping Spectra 



We turn to the case where the spectra of P-space and Q-space overlap, i.e., 
where |-Eo| < 2 A. Here the use of perturbation theory may be somewhat doubt- 
ful but seems justified by the results. We start from the perturbative solu- 
tion (IH]) of the secular equation ([5]). Taking the ensemble average, we change 
the summation over the discrete energies E^j_ into an integration over a con- 
tinuous variable . The negative imaginary increment is needed since the 
spectra overlap, and since the imaginary part of 5a is required to be negative. 
Formally, we obtain the same expressions as in Eqs. flTT]) and f[T^ except that 
I{xq) must be replaced everywhere by 



/-(xo) = /' dx^^^^ . (17) 

J -I X — Xq 



To lowest order in V , this yields 

/A \ r^Eo r^p(xo) 

The real part of the average shift vanishes for Eq = and increases mono- 
tonically in magnitude as Eq moves towards the end points of the semicircle. 
It is negative (positive) for Eq < {Eq > 0), reflecting the effect of the level 
repulsion of the Q-space states below and above Eq. 
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The interaction V mixes the state in P-space with the states in Q-space. Since 
the spectra overlap, that mixing is strong even for small values of the inter- 
action. The degree of mixing is measured by the spreading width which has 
the expected value 27rV^p(-E). For sufficiently small values of T^/X, the prob- 
ability of finding the P-state wave function admixed to the true eigenstates 
of the system is described by a Lorentzian with width [2J. Obviously, the 
spreading width phenomenon has no analogue in the case of non-overlapping 
spectra. The quantity {6a) gives the mean values of shift and width both of 
which actually fluctuate about these mean values. We have not calculated the 
fluctuations. 



4 Generalized Pastur Equation 



In the present Section we use a non-perturbative approach to assess the influ- 
ence of the Q-space states onto the states in P-space. We do so for 1. 
The approach makes use of a generalized form of the Pastur equation. 



4-1 Average Green's Function 



The central element of our analysis is the retarded propagator G{E) (or 
Green's function) of the system, averaged over the GOE. It is deflned as 



Here E is the energy of the system and the plus indicates a positive imagi- 
nary increment. From {G{E)), we find the average level density (or spectral 
function) of the system as 

p=-(l/7r)ImTr(G(P)) , (20) 

where Tr stands for the trace. The function {G{E)) obeys the generalized 
Pastur equation 

{GiE)) = GoiE) + Go{E) (ifGO^(G(E))i/GO^) {G{E)) . (21) 

Here 

GoiE) = i (22) 

is the propagator of the system without any dynamics in Q-space. Eq. (!2T!) 
is obtained by expanding {G{E)) in powers of if^^^, using Wick contraction, 
keeping only leading terms in an asymptotic expansion in inverse powers of 
N (these are the "nested" contributions), and resummation. Alternatively, 
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Eq. (I2T!) can also be derived using supersymmetry [7p] and the saddle-point 
approximation (the saddle-point equation coincides with the generalized Pas- 
tur equation). We use Eq. ([T]) and the definition 



a=^TT{Q{G{E))Q) (23) 



to write Eq. (1211) in the form 

{G{E)) = Go{E) + XaGo{E)Q (GiE)) . (24) 
To solve Eq. fl2^ . we project that equation onto Q-space and obtain 

Taking the trace, we find for a the equation 

A ( QGa(E)Q \ 
"=iV^i l-AaQGo(j;w j 

which can be solved provided QGq{E)Q is known. Inserting the solution into 
Eq. fl2Sl) yields the Q-space projection of the ensemble-averaged propagator. 
The other projections of {G{E)) are easily found, too. We get 



'3<«W>^= l_,.jG„(^)Q QgoWP. (27a) 

P {G(E)) Q = PG„(E)Q(1 + \aQG„(E)Q) , (27b) 
P{G(E))P = PGo(E)P 

+^-PG»(E)Q ^_^^^^^^^^^ QG„iE)P . (27c) 

The various projections of Gq{E) which are needed in the calculation, can be 
obtained using the expansion 



1 °° / 1 \" 



This yields 
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QGo{E)P= — QVP ^ , (29b) 

E+ - QVP ^+ ^ ^ pyg ^+ - ^0 

PGo(i5)g = -;-^—-PVQ ^ , (29c) 

E+ -Ho E+ - PVQ-^QVP 

PGo{E)P = . (29d) 

E+ - PVQ^QVP 

These equations have an obvious physical interpretation. 



4-2 Equation for a 



The Green's function {G{E)) is completely known if we know a, the solution 
of Eq. ( l26l) . To rewrite that equation in a more explicit form we need to work 
out QGo{E)Q. Expanding QGq{E)Q in powers of V , rearranging the series 
and resumming it, we obtain 

QGo{E)Q = -^^^ . (30) 

Here if is a matrix in Q-space given by 

K = Q^ e+^h/ Q ■ (31) 

The matrix K is complex symmetric and can be diagonalized by an energy- 
dependent complex orthogonal transformation. We denote the complex and 
energy-dependent eigenvalues by Kj. The imaginary part of K is negative 
semi definite. Therefore, Im Kj < for all j. The number of nonzero eigenvalues 
Kj is < M. To see this, we use in P-space a basis in which Hq is diagonal and 
has eigenvalues Ej, j = 1, . . . , M. Then, 

K,u = T.V,J^T^V^u. (32) 

This shows that in Q-space and for fixed E, if is a bilinear form in the M vec- 
tors Vfj^j and, therefore, has rank M. We also observe that the imaginary part 
of K is nonzero only if E coincides with one of the eigenvalues Ej. Therefore, 
the M nonzero energy-dependent eigenvalues Kj are, in general, real except 
for a set of discrete points. Using this diagonal form for K, we rewrite Eq. (!26i) 
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as 

N-M A , 1 ^ A 

To discuss Eq. fl33l) we first consider the case where P-space and Q-space are 
uncoupled. Then, K = and Hj = for all j, and Eq. (1551) reduces to the 
well-known saddle-point equation of the GOE which reads 

^ = • (34) 

That equation yields 




and, thus, the semicircle law ([9]) for the normalized average level density of 
the GOE. To compare with Eq. we rewrite Eq. (133|1 as 



a=^A_ + i.f (36) 



- Act (E+ - \a){E+ - Act - 



The saddle-point equation for the GOE is modified. The additional terms re- 
fiect the properties of the Hamiltonian in P-space and its coupling to Q-space. 
While Eq. f lMl) is a quadratic equation in ct and can be solved analytically, 
Eq. fl36l) is of order M + 2 and can only be solved numerically. We use the fol- 
lowing strategy. We first deal with Eq. fl33l) in the case of a one-dimensional 
P-space. Later we show that our main conclusions are not affected as that 
dimension is increased. 



5 One— Dimensional P— Space 



5.1 Basic Equations 



The dimension of Hilbert space is + 1 . As in Section [3l we denote the P- 
space component by the index zero while the indices in Q-space run from 1 
to A^. We write (-ffo)oo = Eq and Vq^ = V^o = We rotate the system in 
Q-space, exploiting the orthogonal invariance of the GOE, to have the vector 
V, initially with components V^, point in the direction of the unit vector in 
the A^-direction. We denote by V the length of that vector. We define 

- E^. ^''^ 
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and have from Eq. (133 



a = X 



N -1 



+ 



A 



1 



N E+ -Xa N E+ -Xa-KQ 
For the projections of Go{E) we obtain 



(3J 



{QGo{E)P) 



fj.1 — Ofj.N- 



-V- 



E+ - Kq E+ - Eq 

{PGo{E)Q),, = 6,M -. ^ - V ^ 



{PGo{E)P) 



E+ - Eq E^ 
1 



Kq 
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E+ -Eq 



E+ 



(39a) 
(39b) 
(39c) 
(39d) 



Eq. (!38l) is the equation we explore numerically in the next Section. We are 
mainly interested in the average level density as given by Eq. (l20l) . To this end 
we need aside from a (which according to Eq. (1231) determines TtQ{G{E))) 
also r = XP{G{E))P. From Eqs. and ([39]), we find 



T = XP {G{E))P 



X{E - Xa) 



{E+ - Eo){E+ - Xa) -V^ 



(40) 



5.2 Perturbative solution of the Pastur equation 



To gain insight into the nature of the solutions, it is instructive to solve the 
generalized Pastur equation for a perturbatively. We write Eq. fl38l) in the 
form 



A 



a 



E-Xa 



N{E-EQ){E-Xa)-V^ ' ^^^^ 
According to Eq. fl20|) . the energies for which the cubic equation fHTj) possesses 
a pair of complex conjugate solutions a (or, equivalently, a single real solution) 
define the spectrum of H. 



5.2.1 Nonoverlapping Spectra 

When Eq lies outside the semicircle, we expect that the spectrum of if consists 
of two disconnected pieces: One piece should more or less coincide with the 
range — 2A < ii^ < 2A of the GOE spectrum. We focus attention on the other 
piece which should cover a small energy interval close to the point Eq. We 
assume Eq < — 2A and \Eq\/X ^ 1. 
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The first factor on tlie riglit-liand side of Eq. (HT!) taken by itself would give 
rise to the unperturbed solution (!35|) which, for Eq < — 2A and \Eq\/X ^ 1, 
is approximately given by X/Eq. We insert this value into the first factor on 
the right-hand side of Eq. (1411) and solve the resulting quadratic equation in 
a. We put E = Eq everywhere except for the term {E — Eq). This yields 



a 



2A 



Ea 



E — Er\ 



± 




E — En 



NEoiE - Eo 



(42) 



The end points of the spectrum are those values of E where the argument of 
the square root vanishes. This yields 




(43) 



E 



To compare this with the result of Section [3l we recall the definition (fTOj) . 
the fact that p(0) = A^/(7rA), and the fact that is the average value of the 
V^^'s while V^^ is their sum. Thus, A^V^ = V^^. We see that the center of the 
interval defined in Eq. (1431) coincides with the shift (IT^ while the length of the 
spectrum is smaller by the factor 2A/|-E'o| than the perturbative result (ITB]) . 



As in SectionOwe briefly address the interaction induced between two states in 
P-space labelled a and b by their interaction with the states in Q-space. Again 
we assume that the two states are degenerate with common energy Eq. A little 
algebra shows that the relevant two-dimensional matrix is Nlajs = {Va\Vp) 
where a and (3 run from a to 6 and where we have used the notation of 
Eq. dS]) (our statement agrees with the arguments in Section [3]). We denote 
the eigenvalues of M^^ by and m^. For simplicity and lack of detailed 
knowledge it is often assumed that the coupling of both states to Q-space is 
the same so that rua = mb- In this case, the Pastur equation takes the same 
form as Eq. fHTj) . with V'^ = rria but with 1/A^ replaced by 2/A^. Making that 
same substitution in Eq. we see that the induced level repulsion is of 
order 1/^/N. This is in keeping with the result of Section [3] if we note that 
rria = mb implies (14 IH) = 0. When that inner product differs from zero, rria 
and rrib must necessarily differ. Expanding the terms in the Pastur equation 
in powers of the difference rria — "^b, it is straightforward to show that the 
leading non-vanishing contributions are of order (rria — m^,)^. That shows that 
significant level repulsion sets in only slowly as |Maf,| increases from zero. We 
take this result as a further justification for studying in Section [6] mainly a 
one-dimensional P-space. 



5.2.2 Overlapping Spectra 

When Eq lies inside the semicircle, we expect the spectrum to consist of a 
single stretch of the energy axis but with a strong enhancement of the average 
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level density for E near Eq. We start from the solution (135|) of the unperturbed 
Pastur equation and write 



a = ao + 6a . (44) 
Inserting this into Eq. (14T]) leads to 



A 

^ N{E- 2Xao - A 6a) [{E - Eq) [E - Xao - X 6a) - V^] ' ^^^^ 

This expression suggests that 6a is of order Hence, we neglect 6a on 

the right-hand side. This approximation does not take into account the shift 
of the end points of the GOE spectrum due to the presence of the state in 
P-space. An improved approximation is introduced below. We find 

_ A 1/2 

" N{E- 2Xao)[{E - Eq){E - Xao) - V^] ' ^^''^ 



(47) 



Likewise we get for r from Eq. fl40|) 

X{E- Xao) 



{E-Eo){E-Xao)-V^ ' 
It is convenient to introduce dimensionless variables. We define 



E Eo 
^=2A' ^°=2A' ^=P- 



As usual we consider the retarded Green's function and, therefore, choose the 
negative sign on the right-hand side of Eq. ( 135|) . We obtain 



6a = -i^ x-iVT^ 

4A^ VT^[x - xf'^' + {i/2)^Vl^] 

and 

11 

T = , . 50 

2 - 7 X - + (i/2)7Vl - ^2 
Both 6a and r display a Breit-Wigner resonance at the shifted energy 

„shift _ _ „ I 7 /p.. N 

Xg — — Xo + -Xo , [Qi) 

Z — 'J z 



with a width 7VI — x"^ where 



We note that 7 is a nonlinear function of 7. To the best of our knowledge, 
that non-linearity has never been taken into account in applications of the 
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doorway-state model to data. The singularity of 7 at 7 = 2 (which corresponds 
to r-*^ = 4A) occurs when the spreading width equals the diameter of the 
semicircle and is, thus, far beyond reasonable applications of the model. To 
lowest order in 7, both shift and width agree with the perturbative result of 
Section [31 



The contribution to the average level density stemming from the P-space state 
is given by the sum of the imaginary parts of N5a and of r. Upon integration 
over all energies these should add up to unity. This is not the case, however, 
because our approximation fails at the end points of the semicircle. Thus, the 
present approximation is useful only if the distance of Eq from the closest end 
point is large compared to V'^/X. 

We improve on this approximation by writing the generalized Pastur equa- 
tion fHTl) in the form 



a{E -Xa) = X 



1 + 



N {E ~ Eo){E - Xa) ~ 



(53) 



We consider this as a quadratic equation in a with known right-hand side. 
This yields 



a 



E 
2X 




+ 



1 



N {E - Eo){E - Xa) - 



(54) 



An approximate solution is obtained by assuming that ^ 1 and substituting 
for a on the right-hand side the unperturbed solution fl5Sl) . Then 



a 



E 
2X 



\ 




+ 



y2 



N(E- Eo){E/2 + iX^l - {E/{2X)y) - ' 

(55) 

For Eq k, 2A near the end point of the semicircle and E > Eq ^ V"^/A, 
the last term under the square root gives a positive contribution: Because 
of level repulsion, the end point of the spectrum is pushed away from the 
center. The argument applies correspondingly when Eq ^ — 2A. The resonance 
contribution can be discussed along similar lines as before. 



6 Numerical Results 



The results of our numerical calculations lead to a deeper understanding of 
the theory developed in the previous Sections. We solve the generalized Pastur 
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N 
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Y=0.5,y„=-2.01 

ky^>-y J =0.25 
perl 



1000 10000 le+05 

N 



Fig. 1. The shift of the single P-space state (M = 1) due to the couphng with 
Q-space states versus N for several values of the coupling strength 7 and for two 
positions yo of the state. 



equation exactly. We use dimensionless variables defined as 



E 



Eq 

T 



7 



A2 



(56) 



We recall that the end points of the semicircle are located at ±2. The defini- 
tion (156|) implies that 7 = r-'^/(2A). Therefore, physically reasonable values of 
7 obey 7 < 1/2 or so, and we have restricted ourselves to that range. 
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Fig. 2. Ratios of exact and perturbative results for the shift for ^ 1 versus 
distance from the centre of the semicircle. 



6.1 Non-overlapping Spectra 



We first consider the case of a one-dimensional P-space hosting the system's 
ground state. The numerical calculation yields the end points of that part of 
the spectrum which is due to the presence of the P-space state. From here 
we calculate the center of the spectrum and its length. We call the difference 
between the unperturbed position and the center of the spectrum the shift and 
the length of the spectrum the fluctuation and compare both values with the 
perturbative estimates of Sections [3] and [5l We must keep in mind, of course, 
that the present definitions of shift and fluctuation differ from the ones used 
in Section [31 

In Fig. [1] the downward shift of the P-space state due to the coupling with the 
chaotic Q-space states is displayed versus N for several values of the cou- 
pling strength 7. For the sake of comparison, the A^-independent perturbative 
result I {no) — yolpert is also given in each panel. In all cases, the downward shift 
increases monotonically with the number N of states in Q-space and reaches 
saturation at about = 10'*. In the sequel, "shift" is used for that asymptotic 
value. 

In the Figure we display two cases: 

a) The P-space state lies far from the Q-space spectrum. The downward shift 
increases monotonically and almost exactly linearly with 7. Indeed, in going 
from 7 = 0.05 to 7 = 0.5 the value grows by an order of magnitude. The 
predictions of perturbation theory are in order and, notably, turn out to be 
always smaller in magnitude than the ones stemming from the generalized 
Pastur equation albeit by a small amount (roughly by about 4 6% for A^ 
not too large). 
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Fig. 3. The fluctuation Ay (square root of the variance) of the position of the single 
P-space state versus for several values of the coupling strength 7 and for two 
positions i/q of the state. We have scaled Ay with N^^"^ . 

b) The P-space state lies just below the semicircle. Again the shift increases 
monotonically with 7 but now perturbation theory badly underestimates 
the shift, especially for small values of 7. 

A summary of our results is given in Fig. [2l For two values of 7, we show the 
ratio of the exact result over the perturbative result for the shift versus distance 
from the centre of the semicircle. We recall that in our units, the radius of the 
semicircle is two. The Figure indicates the goodness of perturbative results. 
For reasonable values of 7 < 1/2, lowest-order perturbation theory furnishes 
reliable results if the distance of the P-state from the edge of the semicircle 
is larger than the radius of the semicircle. 

The fluctuations of the ground-state energy are displayed in Fig. [3l The per- 
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Fig. 4. Ratios of exact and perturbative results for shift and fluctuation. 

turbative result is also given in each panel. In all cases, the fluctuations reach 
saturation at about = 10^. In the sequel, "fluctuation" is used for that 
asymptotic value. 

a) When the P-space state is far from the semicircle we find again that the 
fluctuations grow monotonically with 7, and that perturbation theory is 
reliable but, notably, yields results larger than the ones stemming from the 
Pastur equation (by about 10 -=-15% when 7 is small, up to about 18% when 
7 is large; these estimates depend, of course, upon the definition of Ay). 

b) Perturbation theory fails when the P-space state lies close to the semicircle 
yielding values much smaller than those inferred from the Pastur equation. 

Our results are summarized in Fig. H] where we show the ratios of the exact 
values over the perturbative results for the fluctuations (upper panel) and 
for the shifts (lower panel) versus the coupling strength 7. We observe that 
the discrepancy between perturbative and exact results is decreasing with 
increasing 7. We have no explanation for this observation. 

To assess the validity of our results for realistic cases, we must address the 
case of P-spaces of dimension > 1. We do so for a two-dimensional P-space. 
We assume that the two P-space states are not coupled to each other, and 
that both have the same average coupling to the Q-space states. The results 
for this case are displayed in Fig. [51 We show versus the average shift 
and the fluctuation of the lowest P-space state located originally at i/q = 
—5 for four cases: (i) P-space is one-dimensional (M = 1), (ii) P-space is 
two-dimensional and the second state is originally located at yi = —2.01, 
(iii) P-space is two-dimensional and the second state is originally located at 
yi = —4.9, (iv) P-space is two-dimensional and the second state is originally 
located at 2/1 = —4.99, i.e., in the immediate vicinity of the first state. 
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Fig. 5. Influence of the induced level repulsion between two states in P-space on 
shift and fluctuation versus A^. 

We observe that the downward shift of the lower of the two P-space states 
increases as the distance to the higher P-space state is reduced. In fact, the 
shift now arises from the combined action of the coupling with the Q-space 
and the induced coupling to the other state in P-space. However, this effect 
is visible only for N < 100. For large values of A^, the induced level repulsion 
between the two states in P-space is negligible. This is in accord with the 
results of Sections [3] and [5l The fluctuations also increase as the two states of 
the P-space come closer to each other but again this effect is visible only for 
< 100. 
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6.2 Overlapping Spectra 



We consider a one-dimensional P-space with Eq embedded into the GOE 
spectrum. We display the imaginary parts of the spectral functions Ar and Xa 
as defined in Eqs. (HUj) and In addition we also show these two functions 
for the case where Eq lies far outside the semicircle. 

Figs. [6] and [7] show the imaginary parts for two values of the coupling strength 
7 and for several values of the initial position i/q of the P-space state. In 
each panel we display the results for several choices of the dimension of 
Q-space. We see that in the weak-coupling regime (7 = 0.2) and for i/q = 0, 
the P-space spectral function r is well represented by a Lorentzian. As 
moves towards the end point of the GOE spectrum, the Lorentzian becomes 
increasingly distorted. With the Lorentzian normalized to tt the product of 
height and width (computed at half its maximum) should be 2. For the figure 
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Fig. 6. The imaginary parts of At and of \a for two different coupling strengths 
and for several positions yo of the P-space state. 
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Fig. 7. Same as Fig. [6] but for i/q close to the end point of the semicircle. The four 
panels at the bottom show the values of Ar and of Xa for the case of non-overlapping 
spectra. 
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shown that product is larger than 2 (the more so the more yo approaches the 
center of the semicircle). As a consequence of level repulsion, the imaginary 
part of a is not confined to the interior of the original GOE semicircle. For i/q 
near the center of the semicircle, the dilatation of the semicircle gets smaller, 
however, with increasing A^. The situation changes when i/q is close to the end 
point —2. 

In the strong-coupling regime (7 = 0.5) and for yo near the center of the semi- 
circle, the imaginary part of r resembles a Gaussian more than a Lorentzian. 
Remarkably, the peak height is still given by I/7, irrespective of the strength 
of 7. Again, the function becomes strongly distorted as yo approaches the 
end point. As a consequence of level repulsion, the total spectrum may even 
develop two branches although the P-space state originally lies within the 
semicircle. 

It is of interest to compare the perturbative solution (!55|) of the Pastur equa- 
tion with the exact one. For two values of yo (position of the P-state) and 
several values of A^, this is done in Figs. [Hand O We note that for A^ S> 1, the 
perturbative solution becomes amazingly accurate, even in the strong coupling 
limit. However, when the P-space state is very close to the edge of the semi- 
circle {yo = —1.99), perturbation theory appears to be unable to reproduce 
the position of the peak outside the semicircle. 

We have used the exact results partly displayed in Figs. [6] and [7] to check 
whether the imaginary parts of Na/n and r/vr add up to A^ + 1 as they should. 
We found this indeed to be the case to within about 0.1 percent. This is the 
expected numerical accuracy of our results. Remarkably, this outcome is valid 
no matter what the value of A^: We thus conclude that Pastur equation, while 
missing contributions of the order of 1/A^, yet preserves the normalization of 
the density of states. 



7 Summary and Conclusions 

We have given a comprehensive treatment of the interaction between regular 
and chaotic states. We have modelled the latter in terms of the GOE, the 
Gaussian Orthogonal Ensemble of Random Matrices. Using the limit of large 
matrix dimension for the GOE, we have derived a generalization of the Pastur 
equation for the ensemble average of the Green's function of the system. That 
non-perturbative equation furnishes exact results via numerical calculations. 
Using that approach we have shown that the coupling induced between two 
regular states by their interaction with the chaotic states is very small and, in 
fact, negligible. This is true provided both states are coupled equally strongly 
to the chaotic states. For that reason, we have focused attention throughout 
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Fig. 8. AIm(o" — (To) as obtained from the exact solution of the Pastur equation (soHd 
line) and from the perturbative solution ([55]) (dashed line) for yo = —1, 7 = 0.5 and 
several values for A^. 
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Fig. 9. Same as Fig. El but for yo = —1-99. 
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most of the paper on a single regular state interacting with a large number of 
chaotic states. 

To gain an analytical understanding of the behavior of the system we use 
approximate results. To this end we employ two types of perturbation theory: 
Canonical perturbation theory for the secular equation for the eigenvalues 
of the Hamiltonian, and a perturbative treatment of the generalized Pastur 
equation. 

We have considered two cases: The regular state lies either outside or inside 
the GOE spectrum. In the first case, we calculate the mean value and the 
variance of the shift, or another measure of its spread, as mean values over the 
ensemble. The two perturbative approaches yield essentially identical results 
and agree with the exact solution if the distance between the energy of the 
regular state and the end point of the GOE spectrum is sufficiently large. For 
realistic coupling strengths that means the distance is similar to the radius of 
the semicircle. Due to the interaction with the chaotic states, the spectrum 
of a set of regular levels will get quenched. The quenching is not uniform but 
grows with decreasing distance from the semicircle. 

The interaction not only alters the position (and the wave function) of the 
regular state but also affects the shape of the GOE spectrum. The ensuing 
deformation is particularly relevant when the regular state lies within the GOE 
spectrum and close to one of its end points. The state itself gets shifted, and 
acquires a spreading width. Ganonical perturbation theory, based on a power 
series expansion in the strength of the coupling between regular and chaotic 
states, yields estimates for these quantities. The perturbative solution to the 
Pastur equation does not rely on the smallness of that strength and displays a 
non-linear dependence of the spreading width on the strength. We believe this 
to be an interesting effect which does not seem to have been taken into account 
previously and which deserves further study. The shape of the spectrum due 
to the regular state is Lorentzian for small couplings and gradually changes 
into a Gaussian as the coupling is increased and becomes comparable to the 
radius of the semicircle. The shape change of the chaotic spectrum is not 
accessible to canonical perturbation theory but can be estimated using the 
modified perturbation theory for the generalized Pastur equation. Comparing 
the results of the latter with the exact ones, we find that the perturbative 
solution is excellent for ^ 1. When the energy of the regular state is close 
to the edge of the GOE spectrum, the regular state is pushed outside the GOE 
spectrum as a consequence of level repulsion for sufficiently strong coupling. 
In that case the spectrum develops two branches. The perturbative approach 
is unable to account quantitatively for the part of the spectrum outside the 
semicircle. 

In conclusion: We have derived and solved exactly the generalized Pastur 
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equation for a complex system with strong interactions. We have shown that 
perturbation theory, either canonical or on the Pastur equation, is valid. We 
may view our dynamical system as a generic model for a complex, strongly 
interacting many-body system. In the light of such a view, our results are 
remarkable. They convey the huge simplification of the many-body problem 
induced by the concept of ensemble averaging. 



Appendix 



To calculate the shift of the unperturbed energy £"0 we have to evaluate the 
ensemble averages of the first and of the second term on the right-hand side 
of Eq. ([H]). For the first term we obtain 

I N y2 \ / ^ 1 \ r 1 ^ 

(57) 

The expression I]^i((5(-E' — E^)) is the average level density of the GOE and 
given in Eq. iQ. With the dimensionless variables x = E/2X and xq = Eq/2X 
we get 

where /(xq) is defined in Eq. f|T2l) . As for the second term on the right-hand 
side of Eq. ([8]), we keep only the leading contributions in an expansion in 
inverse powers of N. Thus, in the identity (V^V^) = (V^)^(l + 25^,^) the 
contribution from the Kronecker delta is negligible. Averaging over the V^'s 
yields 

/ N y2 N 1/2 \ 22/^ 1 ^ 1 \ 

v?i Ei^o l^AEu-E,fi = \5 e;:^ S (E^-Eoy) ■ 

(59) 

The two-level correlation function decreases as a Bessel function over distances 
large with respect to the level spacing d, but small with respect to A ~ Nd [7] 
and can be neglected in a calculation keeping only the leading order in 
Thus we get 

For the shift this gives Eq. (|TT1) . The variance of 6a is obtained similarly. 
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